Are networks with more edges easier to synchronize? 
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Abstract. In this paper, the relationship between the network synchronizability and the 
edge distribution of its associated graph is investigated. First, it is shown that adding one 
edge to a cycle definitely decreases the network sychronizability. Then, since sometimes 
the synchronizability can be enhanced by changing the network structure, the question of 
whether the networks with more edges are easier to synchronize is addressed. It is shown by 
examples that the answer is negative. This reveals that generally there are redundant edges in 
a network, which not only make no contributions to synchronization but actually may reduce 
the synchronizability. Moreover, an example shows that the node betweenness centrality is 
not always a good indicator for the network synchronizability. Finally, some more examples 
are presented to illustrate how the network synchronizability varies following the addition 
of edges, where all the examples show that the network synchronizability globally increases 
but locally fluctuates as the number of added edges increases. 
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1 Introduction and problem formulation 

Systems composing of dynamical units are ubiquitous in nature, ranging from physical to 
technological, and to biological fields. These systems can be naturally described by net- 
works with nodes representing the dynamical units and links representing the interactions 

*This work is jointly supported by the National Natural Science Foundation of China under grant 
60674093, the Key Projects of Educational Ministry under grant 107110 and the City University 
of Hong Kong under the Research Enhancement Scheme and SRG grant 7002134. 



among them. The topologies of such networks have been extensively studied and some 
common architectures have been discovered [HE]- The small-world property, for example, 
characterized by short average distance and high clustering among nodes, is one of the most 
common properties shared by many real networks [3]. More significantly, many networks 
show high heterogeneity of node connectivity, which typically possesses a power-law distri- 
bution, named scale-free networks p[]. It is known that these topological characteristics have 
strong influence on the dynamics of the structured systems, such as epidemic spreading, 
traffic congestion, collective synchronization, and so on [5116]. From this viewpoint, system- 
atically understanding the network structural effects on their dynamical processes is of both 
theoretical and practical importance. 

Synchronization behavior, in particular, as a widely observed phenomenon in networked 
systems has received a great deal of attention in the past decades El El [101 El [El ESI [HI 
[T5l [TBI [T71 [18] . Oscillator network models have been commonly used to characterize syn- 
chronization behaviors. In this setting, a synchonizability theorem provided by Pecora and 
Carroll [19] indicates that the collective synchronous behavior of a network is completely de- 
termined by the network structure, independent of individual oscillator dynamics, provided 
that the network coupling strength satisfies some strong conditions. In this framework, it has 
been found that, compared to regular lattices, small- world networks have remarkably better 
synchronizability [20]. In contrast to small- world networks, scale- free networks tend to in- 
hibit synchronization, although they have much shorter average distances than small-world 
networks [21] which are generally deemed to be advantageous for synchronization. Therefore 
the node betweenness centrality was provided as a good indicator to the network synchro- 
nizability [22]. Since the synchronizability is correlated with many topological properties, a 
natural question is which property is the most significant to the synchronizability? Donetti 
et al. [23] gave a good answer to this question by an optimization argument. They pointed 
out that a network with optimized synchronizability should have an extremely homogeneous 
structure, i.e., the distributions of some fundamental topological properties should be very 
narrow. Their work provides a constructive approach to the issue of networked synchroniza- 
tion, making a big progress in this area. However, some issues still remain unclear, e.g., 
what is the most important topological property for the synchronizability? And what is the 
effect of the connectivity density on the synchronizability? In particular, as admitted by the 
authors, this approach cannot theoretically guarantee to find the optimal solution. 

Motivated by the above works, this paper focuses on the relationships between the net- 
work synchronizability and the edge distribution of the associated graph. The effects of the 
connection patterns of graphs on the synchronizability are analyzed both theoretically and 
numerically. It is found that adding an edge to a cycle of size N > 5 definitely decreases 
the network synchronizability, but the synchronizability may be improved by changing the 



cyclic structure. However, a further example shows that, by arbitrarily optimizing the net- 
work structures, networks with more edges are not necessarily easier to synchronize. This 
implies that there are redundant edges in the network with respect to synchronization. 

Consider a dynamical network consisting of N coupled identical nodes, with each node 
being an n-dimensional dynamical system, described by 

N 

Xi = f{xi) - c} j a ij H(x j ), i = 1,2, •• • , TV, (1) 

where Xi = (xn,Xi2, ■ ■ • , Xi n ) G M n is the state vector of node i, /(■) : R n — > W 1 is a smooth 
vector-valued function, constant c > represents the coupling strength, H(-) : M. n — > M. n is 
called the inner linking function, and A = (%-)jVxjv is called the outer coupling matrix or 
topological matrix, which represents the coupling configuration of the entire network. This 
paper only considers the case that the network is diffusively connected, i.e., A is irreducible 
and its entries satisfy 

N 

aa = - aij, % = 1,2, • • • ,N. 

Further, suppose that, if there is an edge between node i and node j, then = = — 1, 
i.e., A is a Laplacian matrix. In this setting, is an eigenvalue of A with multiplicity 1, and 
all the other eigenvalues of A are strictly positive, which are denoted by 

= Ai < A 2 < A 3 < • • • < X N . (2) 

The dynamical network (pQ) is said to achieve (asymptotical) synchronization if x\(t) — > 
X2(t) xpf(t) — > s(t), as t — >• oo, where, because of the diffusive coupling configura- 

tion, the synchronous state s(t) G M. n is a solution of an individual node, i.e., s(t) = f(s(t)). 

It is well known that the eigenratio r(A) = of the network structural matrix A charac- 
terizes the synchronizability. The larger the r(A) is, the better the synchronizability will be. 
The enhancement of the network synchronizability and the relationships between r(A) and 
the network structural characteristics such as average distance, node betweenness, degree 
distribution, clustering coefficient, etc., have been well studied [2H [25j [261 123 [281 [29]. In 
particular, the graph-theoretical method was used to discuss the network synchronizability 
in [2H [2S] ■ This paper further investigates the relationship between the network edges and 
its synchronizability by graph-theoretical tools. 

Throughout this paper, for any given undirected graph G, eigenvalues of G mean eigen- 
values of its corresponding Laplacian matrix. Notations for graphs and their corresponding 
Laplacian matrices are not differentiated, and networks and their corresponding graphs are 
not distinguished, unless otherwise indicated. 



2 Adding one edge to a cycle 



It has been shown by examples that adding edges can either increase or decrease the network 
synchronizability [25], and it was found [30] that in scale-free networks where the nodes are 
coupled symmetrically, if some overloaded edges are removed, the network will become more 
synchronizable. 

In this section, consider adding one edge to a given cycle with N (TV > 4) nodes. In this 
case, the added edge definitely decreases the sychronizability. To show this, the following 
lemmas are needed. 

Lemma 1 [311 132] For any given connected graph G of size N, its nonzero eigenvalues 
indexed as listed in ([2]) grow monotonically with the number of added edges; that is, for any 
added edge e, A;(G + e) > Xi(G), % = 1, • • • , N. 

Lemma 1 shows the eigenvalue changes of graphs due to the addition of edges, but it does 
not show any information about the eigenratio r(A). Therefore, this eigenratio needs to be 
studied in more detail. 

Lemma 2 [271 13T] For any given connected graph G of size N, its largest eigenvalue A at 
satisfies Xn > d max + 1, with equality if and only if d max = N — 1. Further, if G is not a 
complete graph, then the smallest nonzero eigenvalue of G satisfies A2 < d min . Here d max 
and d m i n denote the maximum and minimum degrees of G. 

Lemma 3 [27J For any cycle Cn with N (> 4) nodes, its eigenvalues are given by (jli, ■ ■ ■ , ^n 
(not necessarily ordered as in (T5])) with fxi = and 

sm(f) 

Lemma 4 Given a connected graph G, if the multiplicity of its smallest nonzero eigenvalue 
A 2 is larger than or equal to 2, then adding one edge to G can not change this eigenvalue, 
i.e., X 2 (G + e) = X 2 (G). 

Proof This lemma follows from the fact that rank (A2/ — (G + e)) < rank (A2/ — G) + l. 

□ 

By the above lemmas, one can get the following result for cycles. 

Theorem 1 For any cycle CV with iV > 4 nodes, adding one edge can not enhance its 
synchronizability r(Cjv); specifically, one has r(C^ + e) = r{C^) and r(CV + e) < r(Cjv) 
(N > 5). 

Proof r(C*4 + e) = r{Ci) holds obviously. For the case of > 5, by Lemma 2, one has 
Xn{Cn + e) > 4. But by Lemma 3, Xn{Cn) < 4. And Lemma 3 shows that the multiplicity 
of the smallest nonzero eigenvalue A2 of Cn is 2. By Lemma 4, A2(Ctv + e) = A2(Cat). 
Therefore, r(C N + e) < r(C N ) for all N>h. □ 



Theorem 1 shows that adding one edge a cycle with N > 5 nodes definitely decreases the 
network sychronizability, as shown by the two examples in Figs. 1-5. By simple computation, 
one obtains that r(C 5 ) = ±ff§ = 0.3820 and r(C 5 + e{l,3}) = ^f§ = 0.2993 < r(C 5 ); 
r(C 6 ) = \ = 0.25, r(C 6 + e{l,3}) = T ± m = 0.2265 < r(C 6 ) and r(C 6 + e{l,4}) = | = 0.2 < 
r(C 6 ). 
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Fig. 3 Graph C 6 Fig. 4 Graph C 6 + e{l, 3} Fig. 5 Graph C 6 + e{l, 4} 

From the above two examples, one can find that the synchronizability of cycles strictly 
decreases if only one edge is added, and the results vary depending upon where the edge 
is placed, e.g., r(C 6 + e{l,3}) > r(C 6 + e{l,4}). Considering the optimization of network 
structures, r(C 6 + e{l, 3}) is still not the best one among all graphs with 7 edges connecting 
6 nodes in a cycle, as demonstrated in the next section. 

3 Changing the network structure to enhance its syn- 
chronizability 



It is shown in the above section that adding one edge to a cycle decreases its synchronizability. 
A further question is whether the synchronizability can be enhanced by changing the network 



structure after edge addition. The answer is 'yes' in some cases. For example, one can change 
C5 + e{l,3} to C^ as in Fig. 6, and C6 + e{l,3} to Cq as in Fig. 7. Then, r(Cs ) = | = 0.4 
and r(C 6o ) = ±||f = 0.2684. 




Fig. 6 Graph C 5o Fig. 7 Graph C 6o 

Comparing with the graphs in Figs. 1-5, one can see that both the synchronizabilities 
of C^ and Cq q have been improved. In fact, two cycles share a common edge in Figs. 
2, 4 and 5. In this case, generally the betweenness centrality is large, or the node-to- 
node distances are not homogeneous. In comparison, the network structural characteristics 
are more homogeneous in Figs. 6 and 7. This is consistent with the result of [23]. For 
simple graphs with a few nodes and edges, as those shown above, one can compute their 
eigenvalues to find a good structure for the synchronizability. However, for a general graph, 
how to optimize the network structure toward the best possible synchronizability? Some 
optimal rules are provided based on an optimizing algorithm in [23J: to have homogeneous 
degree, node distance, betweenness, and loop distributions. But these rules are observed 
from simulations, theoretical proofs are not available by now. And, sometimes, these rules 
are contrary to each other. For example, comparing C§ + e{l, 3} with Cq + e{l, 4}, one can 
find that the cycle of C 6 + e{l, 4} is more homogeneous, but the average node distance of C 6 + 
e{l,3} is smaller. It seems that the importance of these rules should be ordered. Although 
some rules are provided in [23], optimizing the network structure for better synchronizability 
is still a hard problem, since it is possible that the optimizing algorithm converges to a 
suboptimal solution. 

Other than the rules for optimization, complementary graphs can be used to characterize 
the network synchronizability [25J. For a given graph G, the complementary graph of G, 
denoted by G c , is the graph containing all the nodes of G and all the edges that are not in 
G. For eigenvalues of graphs and complementary graphs, the following lemma is useful (see 
[3~T| [32] and references therein) . 

Lemma 5 For any given graph G, the following statements hold: 



(i) Ajv(G), the largest eigenvalue of G, satisfies Xn(G) < N. 

(ii) Atv(G) = N if and only if G c is disconnected. 

(iii) If G c is disconnected and has (exactly) q connected components, then the multiplicity 
of \ N (G) = N is q - 1. 

(iv) \ l (G c ) + \ N _ i+2 (G) = N, 2<i<N. 

For example, the complementary graph of C^ is disconnected (see Fig. 8) and the largest 
eigenvalue of C^ is 5, the number of nodes. Its smallest nonzero eigenvalue A2 = 2 can be 
easily obtained by computing the largest eigenvalue of its complementary graph. Further, 
according to the complementary graph, adding one more edge to graph C^ a can not enhance 
its synchronizability. However, if adding two more edges to C 5o , e.g., e{l, 5} and e{3, 5}, then 
the synchronizability increases to r = |. The corresponding complementary graph becomes 
the complementary graph of cycle C4 (see Fig. 9). Cycle C4 and its complementary graph 
are very important in graph theory [32] (see the section below for their further applications). 

For a given graph G, if its complementary graph is disconnected and includes two separated 
graphs (i , and ('•■> . then by Lemma 5 the synchronizability of G is r(G) = v l: ' ;!N ^ A - 1 '• A, "" M( ' :jlf 



iV 



where iV is the number of nodes of G and A max denotes the maximum eigenvalue of the cor- 
responding Laplacian matrix. It is well known that the complementary graphs of bipartite 
graphs are disconnected [25l 133] . so the synchronizability of bipartite graphs can be simply 
analyzed by the above method. Actually, C 5o in Fig. 6 is a bipartite graph. Obviously, better 
understanding and careful manipulation of complementary graphs are useful for enhancing 
the network synchronizability (see the section below for further applications of complemen- 
tary graphs). 




Fig. 8 Graph Cg 



Fig. 9 Graph C% 



4 Are networks with more edges easier to synchronize? 



For a given graph G, let V and £ denote the sets of nodes and edges of G, respectively. A 
graph G\ is called an induced subgraph of G, if the node set Vi of G\ is a subset of V and 
the edges of G\ are all edges in £. In this section, subgraphs and complementary graphs are 
used to discuss network synchronizability. 

In the concern of optimizing network structures, an interesting question is whether net- 
works with more edges are easier to synchronize. In order to answer this question, the 
following lemma is needed. 

Lemma 6 For any given graph G, suppose G\ is its induced subgraph including all nodes 
of G with the maximum degree d max . If G\ includes a cycle Cik with even nodes 2k(k > 2) 
as an induced subgraph, then the largest eigenvalue of G satisfies Xn(G) > d max + 2. 

Proof By Lemma 3, for any cycle Cik with even nodes, its largest eigenvalue is 4. And 
since the degree of every its node is 2, —2 must be an eigenvalue of the adjacency matrix of 
C2k- Let L\ be the sub-matrix of the Laplacian matrix of G associated with all the nodes in 
G\. By the assumption, one has 

(d max + 2)I-L 1 = 2I + A(G 1 )^0, 

where A(G\) is the adjacency matrix of G\. This implies that the largest eigenvalue of G\ 
is larger than or equal to d max + 2. Thus, Lemma 1 leads to the result directly. □ 

Remark 1 Besides Lemma 2, there are few results on lower bounds of the largest eigenvalue 
of Laplacian matrices in graph theory [2H E3j. Since networks with good synchronizability 
always have homogeneous degree distributions, Lemma 6 is very useful for the study of 
network synchronization. 

Theorem 2 For any graph G with 16 edges on 10 nodes, its eigenratio is bounded by 
r(G) < f. 

Proof If the largest node degree of G is d max > 6, then the smallest node degree must 
stisfy d m i n < 2. The conclusion follows directly from Lemma 2. In order to have good 
synchronizability, the degree distribution of G should be homogeneous. Then, first suppose 
that G has 8 nodes with degree 3 and two nodes with degree 4. In this case, by Lemma 2, 
the largest eigenvalue of G is Xio(G) > 5. 

In what follows, consider the largest eigenvalue of the complementary graph G c . By the 
above discussion, G c must have 8 nodes with degree 6 and two nodes with degree 5. Suppose 
G\ is the subgraph of G c composing of 8 degree-6 nodes. By direct computing, G\ must 
have 19 or 20 edges, and the degree of every its node is at least 4. Hence, G\ has 9 or 8 
edges and the degree of every its node is at most 3. If the largest eigenvalue of G\ is 8, i.e., 



G\ is disconnected (Lemma 5), then the largest eigenvalue of G c is larger than or equal to 
8. Therefore, the smallest nonzero eigenvalue of G is A2(G) < 10 — 8 = 2. By the above 
discussion, the theorem obviously holds. Hence, suppose both G% and G\ are connected. 
Then, G\ must have a cycle C4 as an induced subgraph. This holds if and only if G\ has CI 
(see Fig. 9) as an induced subgraph. With only 9 or 8 edges having node degree at most 3, 
drawing G\ directly one can easily reach the conclusion. By Lemma 6, the largest eigenvalue 
of G c must be larger than or equal to 8. Repeating the above discussion concludes the proof. 

When G has 9 nodes with degree 3 and one node with degree 5, the proof can be similarly 
completed. □ 

Remark 2 Theorem 2 shows that there is not a graph G with 16 edges on 10 nodes whose 
synchronizability is r(G) > |. However, there does exist a graph Ti with 15 edges on 10 
nodes whose synchronizability is r(Ti) = | (see Fig. 10), consistent with the result of [23J. 
This clearly shows that networks with more edges are not necessarily easier to synchronize. 
In fact, by the optimal result of [23j, r = | is the optimal synchronizability for graphs 
with 15 edges on 10 nodes. For any graph G with 16 edges on 10 nodes, if both G and G c 
have cycles with even nodes, then by Lemma 6 and Theorem 2, r(G) < ~ = |. Therefore, 
adding one more edge definitely decreases the synchronizability. The existence of cycles with 
even nodes can be easily tested by drawing graphs, so Lemma 6 is very useful for analyzing 
the synchronizability of homogeneous networks. Actually, the graph shown in Fig. 10 is 
quite homogeneous in structure [55] . With one more edge being added, such a structure 
is destroyed. It is therefore easy to understand why adding more edges do not necessarily 
result in better synchronizability. 



1 2 




Remark 3 Fig. 11 shows a new graph T2 with 20 edges on 10 nodes. It also has quite ho- 
mogeneous structural characteristics as discussed in [23] . In fact, the betweenness centrality 



of each node of Ti is 6, larger than that of r 2 , 5. But, the synchronizability of graph r 2 
is worse than that of graph Fx, contrary to the result of [22]. So far, the existing theories 
[2T| [22| 123] [28] can not explain why the synchronizability of I\ is better than that of T 2 . 
This shows the complexity of the relationships between the synchronizability and network 
structural characteristics. Although r 2 has the property of homogeneity, another question is 
whether there exists another graph with 20 edges on 10 nodes having better synchronizabil- 
ity than that of Ti or r 2 ? If the answer is negative, it implies that generally there are many 
redundant edges in a network with respect to its synchronizability. This kind of questions 
are still open today. 
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Fig. 11 Graph T 2 , r(T 2 ) = fgg « 0.382 < r(T x ) = § 



5 Some examples 

In this section, some examples are given to show the changes of the synchronizability versus 
the addition of edges. 

Example 1 The synchronizability changes by adding edges to graphs with cycles are shown 
in Figs. 12 and 13, where their initial graphs are C w and C 50 , respectively, and m a dd denotes 
the number of added edges. The figures in (a)s show the synchronizability changes during 
the process of adding edges with degree homogeneity (i.e., guaranteeing the node degrees 
be as homogeneous as possible during edge-adding). The figures in (b)s show the cases 
corresponding to random edge-adding. Naturally, the corresponding synchronizabilities in 
(a)s are better than those in (b)s, since degree homogeneity is an important property for 
networks to achieve good synchronizabilities. In all graphs, it is shown that the synchroniz- 
ability globally increases but locally fluctuates. According to Theorem 2 and Remark 2, this 
is the expected phenomenon. 




(a) Adding edges with degree homogeneity (b) Randomly adding edges 

Fig. 12 The synchronizability changes of graphs obtained from Cw by adding edges 
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(a) Adding edges with degree homogeneity (b) Randomly adding edges 

Fig. 13 The synchronizability changes of graphs obtained from C 50 by adding edges 

Example 2 The synchronizability changes of graphs obtained from a scale-free graph by 
randomly adding edges are shown in Fig. 14, for which, the same conclusion can be drawn 
as in Example 1. 




Adding edges with degree homogeneity (b) Randomly adding edges 

Fig. 14 The synchronizability changes of graphs obtained from 
a scale-free graph by adding edges 



6 Conclusion 

In this paper, the relationship between the network synchronizability and the edge distribu- 
tion of the associated graphs has been studied. It has been proved that the synchronizability 
definitely decreases if one edge is being added to a cycle with N (N > 5) nodes. However, 
it has also been shown that the synchronizability can be improved by changing the network 
structure. Further, some examples have shown that some networks with more edges, un- 
expectedly, have worse synchronizabilities even if the network structures are in some sense 
optimized. This implies that, for network synchronization, generally there are redundant 
edges, which do not make any contribution to synchronization but may actually destroy the 
synchronizability. In addition, an example of a graph with 20 edges on 10 nodes has been 
provided to show that the existing theories can not explain why it has worse synchronizabil- 
ity than that of a graph with 15 edges on 10 nodes. Some other examples have also been 
given to show that the network synchronizability globally increases but locally fluctuates due 
to edge-adding. According to these results, in practical synchronization problems, the syn- 
chronizability and the number of communication edges should have a coordinative relation. 
And one may utilize the redundant edges to improve robustness or other network properties. 
These kinds of important questions remain open for further research in the future. 
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